Run the Figure4a_4b.do file. More details can be found in the same do file.